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INTRODUCTION 


There have been some studies using 
molecular aspects such as that of Zhang and Shi 
(1993) on phylogeny of rhesus in China and that of 
Melnick et al. (1993) on rhesus from India and 
some of its neighbouring countries. Karanth et al. 
(2008) have researched on molecular phylogeny 
and biogeography of langurs and leaf monkeys 
(Colobine) of South Asia and Zhang and 
Ryder(1998) have worked on old world monkeys 
including Asian Colobines. There have been 
studies on new world monkeys as well (Ascunce 
etal., 2007). Many of these molecular studies have 
been done using blood samples or hair samples 
(Ascunce et al., 2007; Karanth, 2008). Earlier, 
Ascunce et al. (2003) have demonstrated 
amplification of mitochondrial Coll gene from 
DNA extracted from hair samples of some species 
of New World monkeys. However, the present 
study is an attempt on molecular aspects of 
taxonomy on the basis of DNA fingerprinting 
using only hair samples. The objectives of the 
present study included - extraction and 
amplification of DNA from available non-human 
primate hair samples and secondly to process the 
available DNA sequences using appropriate 
softwares to align sequences independently for 
each gene and to conduct phylogenetic and 
molecular evolutionary analysis. 


METHODS 


Initially hair samples obtained from National 
Zoological Collection, Zoological Survey of India, 
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Kolkata were sent to CCMB, Hyderabad for DNA 
fingerprinting. However, as these samples did not 
give good amplifiable DNA, so fresh hair samples 
of monkeys provided by Zoo were sent to CCMB 
for analysis. Out of this also, DNA could not be 
obtained from two samples, so the DNA isolated 
from exhibits A, B, C and F were subjected to PCR 
(Polymerase Chain Reaction)amplification using 
the universal Primers (Verma and Singh, 2003; US 
Patent No. 71,41,364) to generate species specific 
molecular signature. The following protocol was 
used for extraction of DNA from hairs 


Ten to fifteen hair pieces with hair roots about 
0.5 cm were cut and placed into a 1.5 ml eppendorf 
tube. A lysis buffer of 50ul consisting of 1 mm Tris 
pH 8.3, 50 mM KCL and 0.5% Tween was used. 
Besides, 10 1 of 20 g/ml solution of Proteinase K in 
10 mM Tris-HCL (pH 7.5) was also added to the 
above. This was vortexed for 30 seconds. Then, 
ultracentrifuge was done 1t 13000 rpm for one 
second. It was then put in incubator for overnight 
in a 56-60 degree centigrade waterbath. 
Thereafter, it was incubated overnight for 10 
minutes and cooled down to room temperature. 
Again it was ultracentrifuged at 13000 rpm for one 
second. The DNA so extracted was used for PCR. 


Sequences were aligned independently for 
each gene using MUSCLE (Edgar, 2004). 
Phylogenetic and molecular evolutionary 
analyses were conducted using MEGA5 (Tamura 
et al., 2011). Best fit model for nucleotide 
substitution was selected based on minimum 
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Akaike Information Criterion (AIC) value 
(Posada & Crandall, 2001). We constructed 
phylogenetic trees based on maximum 
parsimony (MP) and maximum likelihood (ML). 
Reliability of the phylogenetic tree was estimated 
using bootstrap values run for 1000 iterations. ML 
tree was used for further predicting the molecular 
time scale of species divergence using the 
calibration points 29.2 myr divergence between 
Hylobates and Macaca and 16.7 myr divergence 
between Semnopithecus and Macaca based on time 
tree (Hedges et al., 2006). Hylobates hoolock was 
used as an outgroup. A molecular clock test was 
performed to find out whether the substitution 
rates were uniform (Tamura etal., 2011). 


RESULTS 


The DNA sequences obtained from different 
hair samples are shown in Table 1. Samples A, B 
and F are similar to Macaca mulatta. They are 
monophyletic with a high bootstrap support (Fig. 
1). While sample C belongs to Macaca nemestrina 
which is also supported by high bootstrap value. 
In the cladogram, Hylobates hoolock (an ape) was 
used as an outgroup and other old world monkey 
species were used for comparison and it was 
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found that none of the samples came from 
Semnopithecus entellus (Hanuman langur). 
Maximum likelihood method (Fig. 2) showed 
similar results as Maximum parsimony method 
as in the former method again three samples 
belong to Macaca mulatta and one sample 
belonged to Macaca nemestrina. The sample in 
question was 100% same as Macaca nemestrina 
nemestrina (Gen Bank accession Number 
AF350394). It was first tested whether the 
molecular clock hypothesis is true. Table 2 shows 
that the null hypothesis of equal evolutionary rate 
throughout the tree was accepted (p=0.05). No 
nucleotide divergence was found among samples 
A, B and F as compared to Macaca mulatta 
(Table 3). 


DISCUSSION 


Molecular taxonomy is a useful tool. 
Cytochrome b gene (also used in the present 
study) is widely used in identifying vertebrate 
species. Parson et al. (2000) have used nucleotide 
sequence analysis of cytochrome b gene to 
identify 44 different animal species belonging to5 
major vertebrate groups. Earlier workers have 
used molecular approaches to study primates. 
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Figure 1: Maximum parsimony tree of samples and related species. Sample numbers as per Table 2. There were 310 
bases in final analysis. Bootstrap values are based on 1000 iterations. Hylobates hoolock is used as an 
outgroup. GenBank axccession numbers are provided after the species name. 
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Table - 1: Cytochrome b gene sequence of four samples 


S | Exhibit | BLAST outcome (GenBank Maximum | Sequence 
N Accession No. and species ) identity 
1/A JQ821843 (Macaca mulatta) 100% GGCTTTCATGGGTTATGTTCTCCCATGAGGCCAAATATCATTCTIGGGGAGCAACAGTAATCACAAACCTGC 
TATCAGCAATCCCGTATATCGGAACCAATCTCGTCCAATGAATCTGAGGAGGATACGCCATCGACAGCCCT 
442 ACTCTCACACGATTCTTCACCTTACACTTTATCCTACCCTTCATCATCATCGCCCTCACAACCGTGCACCTA 


CTATTCCTGCACGAAACAGGATCAAACAACCCTTGCGGAATCTCCTCCGACTCAGACAAAATCGCCTTCCA, 
CCCCTACTACACAACCAAAGACATCCTGGGCCTAGTCCTCCTTCTCTTCATCCTAGCAACACTAACACTACT 
CTCACCCAACCTCCTAAACGACCCAGACAACTACATTCCAGCCGACCCATTAAACA 


2|B AY612638 (Macaca mulatta) 100% CACGATTCTTCACCTTACACTTTATCCTACCCTTCATCATCATCGCCCTCACAACCGTGCACCTACTATTCCT 
GCACGAAACAGGATCAAACAACCCTTGCGGAATCTCCTCCGACTCAGACAAAATCGCCTTCCACCCCTAC 
443 TACACAACCAAAGACATCCTGGGCCTAGTCCTCCTTCTCTTCATCCTAGCAACACTAACACTACTCTCACCC 
AACCTCCTAAACGACCCAGACAACTACATTCCAGCCGACCCATTAAACACTCCCCCACATATCAAACCAG 
AGTGATACTTCCTATT 
2e AF35039 100% CACCCTACACTTCATCCTACCCTTCATCATCGTCGCCCTCACAACCGTACACCTACTTITTCCTACACGAAAC 
AGGATCAAACAACCCCTGCGGAATITCCTCCGACTCGGATAAAATCACCTTTCACCCCTATTATACAATCA 
444 (Macaca nemestrina nemestrina) AAGACATCCTAGGCCTAATCCTCCTTCTCTTTGCCCTAACAATACTAACACTATTCTTACCCAACCTCCTAA 


ACGACCCAGACAACTACATTCCAGCTGACCCACTAAATACCCCCCCACATATCAAGCCAGAGTGATACTT 
CCTATTTGCATACACAATCCTACGAT 


4|F JQ821843 (Macaca mulatta) 99% TCTCACACGATTCTTCACCTTACACTTTATCCTACCCTTCATCATCATCGCCCTCACAACCGTGCACCTACTA 
TICCTGCACGAAACAGGATCAAACAACCCTTGCGGAATCTCCTCCGACTCAGACAAAATCGCCTTCCACCC 
447 CTACTACACAATCAAAGACATCCTGGGCCTAGTCCTCCTTCTCTTCATCCTAGCAACACTAACACTACTCTC 


ACCCAACCTCCTAAACGACCCAGACAACTACATTCCAGCCGACCCATTAAACACTCCCCCACATATCAAA 
CCAGAGTGATACTTCCTATTTGCATACACAATCCTACGATCCATCCCCAACAAACTGGGAGGCGTACTAGC 
ACTCTTTCTATCGATCCTCATTCTAGCAGCCATCCCT 
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Figure 2: Molecular time clock for the evolution of Macaca and related species based on maximum likelihood 
method with Tamura and Nai (1993) model of nucleitide substitution and invariant sites. Hylobates hoolock 
is used as an outgroup. GenBank axccession numbers are provided after the species name. Calibration 
points are 29.2 myr divergence between Hylobates and Macaca and 16.7 myr divergence between 
Semnopithecus and Macaca based on time tree (Hedges et al., 2006). 
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Table-2: Results from a test of molecular clocks using the Maximum Likelihood method. The null 


hypothesis of equal evolutionary rate throughout the tree was accepted 
(P =0.05). 


Parameters 


With Clock -1273.248 23 


Without Clock -1259.835 


Table-3: Nucleotide Sequence of Various Macaca Species. 


10 20 30 40 50 60 

jaaa kaas losa. kesan, ioes [Sees [eee [pase [as ana baas! boaa iea 
Macaca_mulatta_JQ821843 CACCUUACACUUUAUCCUACCCUUCAUCAUCAUCGCCCUCACAACCGUGCACCUACUAUUCC 
Macaca_mulat ta_AY612638 CACCUUACACUUUAUCCUACCCUUICAUCAUCAUCGCCCUCACAACCGUGCACCUACUAUUCC 
X442 CACCUUACACUUUAUCCUACCCUUCAUCAUCAUCGCCCUCACAACCGUGCACCUACUAUUCC 
X443 CACCUUACACUUUAUCCUACCCUUCAUCAUCAUCGCCCLUICACAACCGUGCACCUACUAUUCC 
X447 CACCUUACACUUUAUCCUACCCUUCAUCAUCAUCGCCCUCACAACCGUGCACCUACUAUUCC 
Macaca_nemestrina_nemestrina_AF350394 CACCCUACACUUCAUCCUACCCUUCAUCAUCGUCGCCCUCACAACCGUACACCUACUUUUCC 
X444 CACCCUACACUUCAUCCUACCCUUCAUCAUCGUCGCCCUCACAACCGUACACCUACUUUUCC 


Macaca_nemestrina_nemestrina_AF350397 CACCCUACACUUUAUCCUACCCUUCAUCAUCGUCGCCCUCACAACCGUACACCUACUCUUCC 
Macaca_nemestrina_leonina_AF350399 CACCCUGCACUUCAUCCUACCCUUICAUCAUCGUCGCCCUCACAACCGUACACCUACUUUUCC 
Macaca_nemestrina_leonina_AF350398 CACCCUGCACUUCAUCCUACCCULICAUCAUCGUCGCCCUCACAACCGUACACCUACUUUUCC 


Macaca_fascicularis_FJ906803 CACCCUGCACUUUAUCCUACCCUUCAUCAUCAUCGCCCUCGCAGUCGUACACUUACUAUUCC 
70 80 90 100 110 120 
Macaca_mulatta_JQ821843 UGCACGAAACAGGAUCAAACAACCCUUGCGGAAUCUCCUCCGACUCAGACAAAAUCGCCUUC 
Macaca_mulatta_AY612638 UGCACGAAACAGGAUCAAACAACCCUUGCGGAAUCUCCUCCGACUCAGACAAAAUCGCCUUC 
X442 UGCACGAAACAGGAUCAAACAACCCUUGCGGAAUCUCCUCCGACUCAGACAAAAUCGCCUUC 
X443 UGCACGAAACAGGAUCAAACAACCCUUGCGGAAUCUCCUCCGACUCAGACAAAAUCGCCUUC 
X447 UGCACGAAACAGGAUCAAACAACCCUUGCGGAAUCUCCUCCGACUCAGACAAAAUCGCCUUC 
Macaca_nemestrina_nemestrina_AF350394 UACACGAAACAGGAUCAAACAACCCCUGCGGAAUUUCCUCCGACUCGGAUAAAAUCACCUUU 
X444 UACACGAAACAGGAUCAAACAACCCCUGCGGAAUUUCCUCCGACUCGGAUAAAAUCACCUUJ 
Macaca_nemestrina_nemestrina_AF350397 UACACGAAACAGGAUCAAAUAACCCCUGUGGAAUUUCCIUCCGACUCAGACAAAAUCGCCUUJ 
Macaca_nemestrina_leonina_AF350399 UACACGAAACAGGAUCAAACAACCCCUGCGGAAUUUCCUCCGACUCGGAUAAAAUCACCUUU 
Macaca_nemestrina_leonina_AF350398 UACACGGAACAGGAUCAAACAACCCCUGCGGAAUUUCCUCCGACUCGGAUAAAAUCACCUUL 
Macaca_fascicularis_FJ906803 UACACGAAACAGGAUCAAACAACCCUUGUGGAAUUU CCUCCGACUCGGACAAAAUCGCCUUC 
130 140 150 160 170 180 

|--==|---=|-=--[-=-=|---- | ----|---- | ----|---- |---- |----|---- - 
Macaca_mulatta_JQ821843 CACCCCUACUACACAACCAAAGACAUCCUGGGCCUAGUCCUCCUUCUCUUCAUCCUAGCAAC 
Macaca_mulatta_AY612638 CACCCCUACUACACAACCAAAGACAUCCUGGGCCUAGUCCUCCUUCUCUUCAUCCUAGCAAC 
X442 CACCCCUACUACACAACCAAAGACAUCCUGGGCCUAGUCCUCCUUCUCUUCAUCCUAGCAAC 
X443 CACCCCUACUACACAACCAAAGACAUCCUGGGCCUAGUCCUCCUUCUCUUCAUCCUAGCAAC 
X447 CACCCCUACUACACAAUCAAAGACAUCCUGGGCCUAGUCCUCCUUCUCUUCAUCCUAGCAAC 
Macaca_nemestrina_nemestrina_AF350394 CACCCCUAUUAUACAAUCAAAGACAUCCUAGGCCUAAUCCUCCUUCUCUUUGCCCUAACAAU 
X444 CACCCCUAUUAUACAAUCAAAGACAUCCUAGGCCUAAUCCUCCUUCUCUUUGCCCUAACAAU 


Macaca_nemestrina_nemestrina_AF350397 CAUCCCUACUAUACAGUCAAAGACGUCCUAGGUCUAGUUCUCCUUCUAUUCGUUCUGACAAC 
Macaca_nemestrina_leonina_AF350399 CACCCCUACUACACAAUCAAAGACAUCCUAGGCCUAAUCCUCCUUCUCUULUGCCCUAGCAAC 
Macaca_nemestrina_leonina_AF350398 CACCCCUACUACACAAUCAAAGACAUCCUAGGCCUAAUCCUCCUUCUCUULIGCCCUAGCAAC 


Macaca_fascicularis_FJ906803 CACCCCUACUACACAAUCAAAGACAUCCUAGGUCUGAUCCUCCUUCUCUUCAUCCUAGCAAC 
190 200 210 220 230 240 

aes | cman [lee | eset | ena |==== | --== | -—=== | =-=-= |---- | =--= |---- | === 
Macaca_mulatta_JQ821843 ACUAACACUACUCUCACCCAACCUCCUAAACGACCCAGACAACUACAUUCCAGCCGACCCAU 
Macaca_mulatta_AY612638 ACUAACACUACUCLICACCCAACCUCCUAAACGACCCAGACAACUACAUUCCAGCCGACCCAU 
X442 ACUAACACUACUCUCACCCAACCUCCUAAACGACCCAGACAACUACAUUCCAGCCGACCCAU 
X443 ACUAACACUACUCUCACCCAACCUCCUAAACGACCCAGACAACUACAUUCCAGCCGACCCAU 
X447 ACUAACACUACUCUCACCCAACCUCCUAAACGACCCAGACAACUACAUUCCAGCCGACCCAU 
Macaca_nemestrina_nemestrina_AF350394 ACUAACACUAUUCUUACCCAACCUCCUAAACGACCCAGACAACUACAUUCCAGCUGACCCAC 
X444 ACUAACACUAUUCUUACCCAACCUCCUAAACGACCCAGACAACUACAUUCCAGCUGACCCAC 


Macaca_nemestrina_nemestrina_AF350397 ACUAACACUAUUCUCACCCAACCUCCUAAAUGACCCAAAUAACUACAUUCCAGCOGACCCAC 
Macaca_nemestrina_leonina_AF350399 ACUAACACUAUUCUUACCCAAUCUCCUAAACGACCCAGACAACUAUAUUCCAG ‘ACCCAC 
Macaca_nemestrina_leonina_AF350398 ACUAACACUAUUCUUACCCAAUCUCCUAAACGACCCAGACAACUAUAUUICCAGCUGACCCAC 
Macaca_fascicularis_FJ906803 ACUAACACUACUCUCACCCAACCUCCUAAACGACCCAGAUAACUACACUCCAGCCGACCCAC 


| t iM | 


250 

ap amne 
Macaca_mulatta_JQ821843 UAAACA 
Macaca_mulatta_AY612638 UAAACA 
X442 UAAACA 
X443 UAAACA 
X447 UAAACA 
Macaca_nemestrina_nemestrina_AF350394 \JAAAUA 
X444 UAAAUA 
Macaca_nemestrina_nemestrina_AF350397 UAAAUA 
Macaca_nemestrina_leonina_AF350399 UAAAUA 
Macaca_nemestrina_leonina_AF350398 UAAAUA 
Macaca_fascicularis_FJ906803 UAAACA 


! 
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Zhang and Shi (1993) classified rhesus in China 
into 6 sub-species using mitochondrial DNA. 
Work of Karanth et al. (2008) using two unlinked 
nDNA markers (lysozyme and protamine P1) and 
one mtDNA (cytochrome b) marker on langurs 
and leaf monkeys from South Asia support the 
classification of langurs of Indian sub-continent 
(Hanuman, Nilgiri and Purple-faced langurs) in 
the genus Semopithecus and leaf monkeys from 
South-East Asiain the genus Trachypithecus. 


Macaca mulatta and Macaca nemestrina 
diverged around 9 m years ago i.e. mean 
divergence is 9 m years ago (confidence interval 
4-14 m years ago, because one gene only was used 
so confidence interval was very large. 


Ascunce et al. (2007) have used blood and hair 
samples for their molecular study on New World 
Monkey species Alouatta caraya and Karanth 
et al. (2008) have used blood, tissue and hair 
samples. The present work is based exclusively on 
DNA extracted from hair samples. So, the 
approach used in the present study is non- 
invasive based on hair samples which may be 
useful in identification of seized material. Besides, 
it may also be useful in studies involving free- 
living animals (removing the problem of catching 
the individuals for taking out blood samples). 


SUMMARY 


The present work uses a non-invasive 
approach for study of taxonomy of some of the 


Indian non-human primates based on DNA 
sequences obtained through cytochrome b gene 
from hairs. The DNA isolated was subjected to 
PCR (Polymerase Chain Reaction) amplification 
using the universal Primers to generate species 
specific molecular signature. The DNA 
sequences sent by CCMB, Hyderabad were 
analysed . 


Reliability of phylogenetic tree was estimated 
using bootstrap values run for 1000 iterations. A 
molecular clock test to check whether substitution 
rates were uniform showed that the null 
hypothesis of equal evolutionary rate throughout 
the tree was accepted (p=0.05). Maximum 
likelihood method showed results similar to that 
of Maximum Parsimony method. Nucleotide 
substitution sites in cytochorme-b gene among 
different study species were also found out. No 
nucleotide divergence was found among Macaca 
mulatta samples in question. Future studies 
incorporating more genes and larger number of 
samples and species may provide further insight. 


DNA based taxonomy is useful in wildlife 
forensics and non-invasive method is particularly 
useful in free-living and zoo primates as it does 
notinvolve catching them. 
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